Haze of surface random systems: An approximate analytic approach 
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G^ ■ 

■ Approximate analytic expressions for haze (and gloss) of Gaussian randomly rough surfaces for 

' various types of correlation functions are derived within phase-perturbation theory. The approxima- 

, tions depend on the angle of incidence, polarization of the incident light, the surface roughness, a, 

and the average of the power spectrum taken over a small angular interval about the specular direc- 
^ ' tion. In particular it is demonstrated that haze(gloss) increase(decrease) with a/X as exp{—A{a/ X)^) 

and decreases(increase) with a/X, where a is the correlation length of the surface roughness, in a way 
that depends on the specific form of the correlation function being considered. These approxima- 
fT^ ' tions are compared to what can be obtained from a rigorous Monte Carlo simulation approach, and 

I good agreement is found over large regions of parameter space. Some experimental results for the 

angular distribution of the transmitted light through polymer films, and their haze, are presented 
and compared to the analytic approximations derived in this paper. A satisfactory agreement is 
^ ' found. In the literature haze of blown polyethylene films has been related to surface roughness. Few 

authors have quantified the roughness and other have pointed to the difficulty in finding the correct 
roughness measure. 
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PACS numbers: 42.25.-p; 41.20.-q 



I. INTRODUCTION 



Optical properties of polyethylene films have attracted considerable attention due to the importance in applications 
such as packaging. Studies of haze and gloss of films have been carried out addressing the effect of polymer structure [J, 
0, d, 0] , Theological properties 0, Q , additives [1, 0, 11] and processing conditions [1, H, [13, [HI [12, [13] ■ 
' The notion of haze ("cloudiness") is supposed to quantify the ratio between the diffusely reflected or transmitted 
: light to the total light reflected (reflectance) or transmitted (transmittance). To this end, haze of a fllm is defined 
CO ■ as the fraction of transmitted light that deviates from the directly transmitted beam by more than given amount 
' {e.g. 2.5°) [3, [H, [T^]. A similar definition applies for reflection. For thin fllms, haze is recognized to be caused 
mainly by scattering from surface irregularities, in contrast to bulk randomness. Previously, only a few groups have 
^\ reported on the dependence of haze on surface roughness [1, [l^ . [l7j . In a study of different polyethylene materials, 
the bulk contribution to haze of 40/im thick fllms was found to vary in the range 10-30% of the total value [|. 
Two main mechanisms for surface roughness have been identifled [l^: (i) flow- induced irregularities originating from 
the die (extrusion haze), and (ii) protruding crystalline structure such as lamellae, stacks of lamellae or spherulites 
(crystallization haze) . 

Surface roughness is often characterized (in the engineering literature) by the root-mean-square deviation from the 
mean surface height, a. Different roughness generating mechanisms, such as die irregularities and inhomogeneous 
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distribution of additives, will influence the roughness at different length scales. Large-scale trends (compared to the 
wavelength) will not affect the diffuse light scattering and must be eliminated from the analyzes. Implicit difficulties 
in deriving a relevant measure of surface roughness may obscure its correlation with haze. Various techniques have 
been applied to characterize roughness. In a number of recent studies atomic force microscopy (AFM) has been 
applied [1, [13, 0, HO, IMI ■ Many authors, however, associate only qualitative differences in roughness, as observed by 
e.g. AFM, with haze values. Robust methods still have to be developed to extract the relevant roughness measure 
based on a sufficiently high statistics. A discussion of surfaces characterizing and relation to haze is given in 

Sukhadia et al. |22| relate both crystallization haze and extrusion haze to the elastic properties measured as the 
recoverable shear strain of the polymer. For a low level of melt elasticity orientations from the die relax quickly and 
crystalline aggregates form on or close to the film surface. For highly elastic materials the surface roughness was 
attributed to melt flow effects generated at the die exit. Plotting haze versus the recoverable shear strain J2^, based 
on a large data set of blown and cast films, results in a parabolic curve with lowest haze for materials of intermediate 
level of elasticity. 

Studies of electromagnetic wave scattering have a long history, and the effect of surface roughness on the scattering 
has also been studied for many years in optics [H, [13, HI, [13, [H, [H, [13, [U, [H, [H, [13, l35| and more recently, for 
x-ray scattering [H, [13, [H, [la \^ ■ Today it is fair to say that the main features are rather well understood, at least 
for one-dimensional roughness, where one recently has started to address inverse (optical) problems [U, [43|. In the 
case of two-dimensional roughness there are still open questions to be answered, in particular for optical frequencies 
where the dielectric contrast (and therefore the scattering) is the most pronounced. This paper concentrates on the 
range of optical frequencies, even if the theoretical results derived herein hold for any frequency for which the adopted 
approximations are expected to hold. This choice is made due to the fact that the concept of haze (and gloss) mainly 
is used within optics, to the best of our knowledge. 

Different approximations are available for scattering of electromagnetic waves scattered by rough surfaces or bulk 
inhomogeneities, depending e.g . on the typ ical range of roughness/ wavelength ratios (cr/A), scattering angles, and 
electrical conductivity [13, [Ifl [13, [H, l44| . Vectorial or scalar formulations are used, depending on whether or 
not polarization is taken into account. Vectorial formulations, such as the Rayleigh-Rice perturbation theory, are 
considered to be more accurate for smooth surfaces (typically a <^ X) with finite conductivity, and large scattering 
angles. However, the scalar Kirchhoff theory is claimed to be more accurate than Raylcigh-Rice for rougher surfaces, 
and it is easier to handle mathematically ^43;]. 

Some calculations of gloss vs. surface roughness parameters have been published over the last years. Alexander-Katz 
and Barrera [1^ calculated the angular distribution of reflected light, using the scalar Kirchhoff approximation. With 
a Gaussian height distribution and various height correlation functions, the gloss depended on two parameters. It 
increased with decreasing a and increasing correlation length. The sensitivity of gloss to correlation length depended 
on a. In particular, the sensitivity was low for very low and high cr/A values. Wang et al. [20] performed similar 
calculations. With surface parameters (from AFM) and refractive indices as input, the calculations agreed fairly well 
with gloss values in the range 30-70%, measured on PE films. 

Haze has been calculated as function of bulk inhomogeneities [13], but theoretical literature on haze vs. surface 
roughness is sparse, although there are related studies e.g. on radio transmission Willmouth [13] referred to 

unpublished geometrical optics calculations relating film clarity to surface roughness on a scale greater than the 
wavelength of light. In most cases, however, the roughness is on a finer scale, and the calculations must be based on 
physical optics. The scalar Kirchhoff approximation can, in principle, be applied to transmission calculations [13], 
but this is more difficult than the reflection case, due to contributions from two surfaces, possible bulk effects, and a 
high probability of multiple scattering [4^. 

More recently Wang et al. [ll| presented a method for calculating haze. Since the (real) surfaces were spherulitic, 
the combined scattering of the two surfaces could be modelled by Mie scattering theory (valid for a single sphere 
of any size and refractive index). Calculated haze values agreed with experimental data for films of six different 
materials. Furthermore, the model predicted a maximum in haze for a spherulite diameter of 800 nm (A = 550 nm), 
while the clarity decreased monotonously with increasing spherulite diameter. This critical diameter, corresponding 
to maximum haze, decreased with increasing refractive index, but it was insensitive to the volume fraction of spheres, 
in the range studied. Although this model provides insight into the relationship between surface roughness and haze, 
it has some limitations. In particular, it can not account for non-spherical protrusions and lateral correlation. 

In this paper we study the effect of surface roughness on haze by addressing both the height distribution as well 
as the height-height lateral correlations. A motivation has been that problems in deriving experimental measures 
of surface roughness have obscured the correlation with haze [3]- Our main goal is to establish the relationship 
between surface structure — both roughness amplitude and surface lateral correlations — and haze. An analytical 
approximation to haze for Gaussian randomly rough surfaces is derived which is compared with a rigorous Monte 
Carlo simulation approach. Models are compared with experimental results on haze as well as with the angular 
distribution of the transmitted light through films. 



This paper is organized as follows: The following section introduces the scattering geometry to be discussed in this 
work as well as the properties of the surface roughness. In Sec. IIIII the scattering theory and notation to be used in 
the following discussion is presented. The main results of this work starts in Sec. llVl where the definition and analytic 
approximate expressions for haze are presented. This approximation is compared to rigorous computer simulations 
in Sec. |Vl Finally the conclusions that can be drawn from the present work are presented in Sec. IVII 



II. THE SCATTERING GEOMETRY 



The scattering geometry that we will consider in this study is depicted in Fig. [T] In the region z > (^{x) it consists 
of vacuum {eo{uj) = 1) and for z < ({x) of a dielectric characterized by an isotropic, frequency-dependent, dielectric 
function ei{uj). Here ({x) denotes the surface profile function. It is assumed to be a single- valued function of x that 
is differential as many times as is necessary. Furthermore, it constitutes a zero mean, stationary, Gaussian random 
process that is defined by 

iCix)) = 0, (la) 
{C{x)ax')) = a'W{\x~x'\). (lb) 

Here W^(|x|) denotes the (normalized) auto-, or height-height correlation function and will be specified later, a is the 
root-mean-square of the surface roughness, and (•) denotes the average over an ensemble of realizations of the surface 
roughness. For the later discussion one will also need the power spectrum of the surface roughness, defined as the 
Fourier transform of the correlation function, i.e. 

g{\k\) = / dxe-''^Wi\x\). (2) 



In this paper, we will mainly deal with correlation function of the exponential type, W{x) = exp(— |a;|/a), where a is 
the so-called correlation length. For such a correlation function the power spectrum becomes 

The incident wave will be assumed to be either p- or s-polarized, as indicated by the subscript v on field quantities, 
and the plane of incidence will be the xz-plane. Furthermore, the angle of incidence, scattering, and transmission, 9q, 
9s 1 and 9t respectively, are measured positive according to the convention indicated in Fig.[TJ 

In order to demonstrate that the assumption made above for the statistics of the surface roughness is not unre- 
alistic, we in Fig. ^a.) present an AFM measurement of the surface topography of a polyethylene film surface. The 
corresponding height distribution and height-height correlation function that can be obtained on the basis of such 
topography measurements are can be found in Figs. [Hb) and (c). It is observed from these figures that the mea- 
sured surface roughness, to a good approximation, is a Gaussian random process. For this particular example, an 
exponential correlation function is a reasonable, but not perfect, choice for the correlation function. 



III. SCATTERING THEORY 



In this section, elements of the scattering theory that will be useful for the discussion that will follow will be 
presented. A more complete and detailed presentation can be found e.g. in Rcf. 31]. 



A. The reflection and transmission amplitudes 



Due to the one-dimensional character of the surface topography, z ~ C{x), a generic scalar field may be introduced 
that fully can describe, together with Maxwell equations, the electromagnetic field. Such generic field is defined as 

where z/ is a polarization index, and Hy and Ey are used to denote the second component of the electric and magnetic 
fields. The advantage of using the field variable (l>^{x,z,uj) is that the Maxwell (vector) equations satisfied by the 



electromagnetic field are equivalent to the scalar wave equation for (j)^(x, z,lu). To be able to use a scalar equation 
instead of vector equations represents a great simplification of the problem. However, it should be noted that this 
simplification comes about because of C(a;) being one-dimensional, and it doesn't hold true in general. 
For the scattering system considered in this paper, the field can be written as 

when z > laax^^x) and a plane incident wave is assumed, and 

M^,^,^) = ^r,(g|A;)e''™(«)^ (5b) 

when z < mm({z). In writing Eqs. (O, we have introduced R^{q\k) and T^{q\k) — the reflection and transmission 
amplitudes, the lateral momentum variable 

q = ^/e — sin9, (6) 
c 

where e is the dielectric constant of the medium considered, and 6 denotes the angle of incidence, reflection or 
transmission depending on context. Furthermore, in Eqs. ([5]), we have also defined 



"™('?) = { ^ I , (7) 

i^Jq^-Em^, l^l > 

where m — Q corresponds to the medium above the rough surface, and to = 1 to the dielectric medium below. Notice 
that when \q\ < y^e^{uj/c) (non-radiative region), it follows that am{q) = ^em('-^/c) cos 9. 



B. Mean Differential Reflection and Transmission Coefficients 



The mean differential reflection and transmission coefficients, abbreviated DRC and DTC, are two experimentally 
and theoretically accessible quantities frequently used to study the angular distribution of the reflected or transmitted 
light. We will here denote them by {dRu/dOg) and {dT^/dOt) respectively, where subscript v, as before mentioned, is 
a polarization index. The mean DRC is defined as the fraction of the incident power that is scattered by the rough 
surface into an angular interval ddg about the scattering angle Og. The power (energy flux) crossing a plane parallel 
to the xy-plane can be calculated from 

' dxj dy ReiSs)^, (8) 

where S = 1/2 E x H* is the (complex) Poyntings vector [HI for the electromagnetic field (E, H), and (•)j denotes a 
time-average. By using Eqs. ([5]) one finds that the incident power is given by 

Pmc = — ao(fc), (9a) 

Z UJ 

while the scattered power becomes 

Psc = eo^- / ^ao{k)\RMk)\' = / d9sPsc{es). (9b) 

This latter relation implicitly defines the angular dependent scattered power Psdds)- Thus from the definition of the 
differential reflection coefficient it is realized that dR,y/d9s = Psci^s)/ Pine- However, since the surface is randomly 
rough, it is the mean differential reflection coefficient that should be of interest. Such a quantity will be given by 

where we have used (• • • ) to denote the average over surface realizations. In the same way, the mean differential 
transmission coefficient may be defined as {dT^/dOt) = {ptri^t) / Pine) , where ptr{&t) denotes the power transmitted 



into an angular interval about the angle of transmission Ot- Instead of giving the expression for this quantity explicitly, 
we will instead for later convenience, introduce a generic notation for both the mean differential reflection and 
transmission coefficient. This generic quantity will be denoted by {dU/d6) so that 

dU\ [(w"/' in reflection, 

(11a) 

in transmission 

where 9 stands for 0^ and 6t in reflection and transmission, respectively. In general it may be written in the following 
form 




dU \ 1 £m w cos^ ( 



\U{q\k)\'). (lib) 



80 I Li ^/sq I-kc cos 6*0 

In Eq. (jllb[) {/(gj/c) denotes the reflection or transmission amplitudes depending on context, and m takes on the 
values m = in reflection and m = 1 in transmission. 

In theoretical studies it is customary to separate (dU j 80) into two terms — one coherent and one incoherent term. 
That this is possible can be realized from the following (trivial) rewriting 



(\U{q\kf) = \{U{q\k))\' + [(\U{q\kt) \{U{q\k))\^] . (12) 



If this expression is substituted back into Eq. (|llb[) . the first term will give rise to the coherent or specular contribution 
to {811/89), while the last term, within the square brackets, will result in the incohrent or diffuse contribution to 
the same quantity. We will use subscripts coh and incoh, respectively, to indicate these whenever needed. This 
identification follows from observing that in the average {U{q\k)) one only gets contributions from those portions of 
U{q\k) that are in phase from one surface realization to another. 

From Eq. (|llb|) a quantity that defines the fraction of the incident energy that is either reflected or transmitted can 
be defined as 



7r/2 



89 



Notice, that when no absorption takes place in neither media involved, i.e. ImSm — 0, one should have that the 
sum of this quantity in reflection and transmission should add up to one, Lis + hit — 1- This is a direct consequence 
of energy conservation. In practice, however, there will always be some absorption, but for many dielectric media at 
optical frequencies it is a reasonable approximation to neglect it. 



IV. HAZE 



In the optical industry two quantities — haze and gloss — are often used to quantify the visual appearance of 
materials 43]. Gloss, crudely speaking, is related to the amount of light being reflected (or transmitted) into angles 
around the specular direction. This quantity has previously been studied experimentally [5Q] for transparent plastic 
materials, and recently also studied theoretically [4^| . In this latter study, the authors investigated how gloss depends 
on the level of roughness and surface correlations. It was found that the incoherent (diffuse) contribution to the 
scattered light could contribute significantly to the gloss. Haze, on the other hand, measures the fraction of reflected 
(transmitted) light that is reflected (transmitted) away from the specular direction [l^ [l^ . In the former case one 
talks of haze in reflection and in the latter of haze in transmission. Notice, that haze can almost be considered as a 
complementary quantity to gloss. If the haze of a transparent plastic film, say, is large, then an object viewed through 
the film will look unsharp or blurry. It is this kind of visual effect that the haze value is supposed to quantify. 

Even though gloss is important in many applications, we will in the present study concentrate on haze since this 
quantity has not been studied that extensively in the literature by theoretically means. 



A. Definition of haze 



For the purpose of this theoretical study, we will focus on a semi-infinite medium, instead of a film geometry. The 
reason for this choice is purely practical. Let us start by assuming that the incident light is impinged onto the planar 
mean surface at an angle measured counter clock- vise from the normal to the mean surface (cf. Fig.[T|). In case of 



no surface roughness the hght wiU be scattered (or transmitted) in accordance with Snell's law. Hence all the energy 
will be propagate in the directions defined by the angle 



e = arcsin I -^sin6io ) , (14) 



where m = should be used in reflection (for which <d = Bq), and to = 1 in transmission. 

Formally, haze, Tl{6o), is defined [l^, [la, IBH as the fraction of the reflected (transmitted) light that is reflected 
(transmitted) into angles lying outside the angular interval (9^,9+) with 

0±=e±A9, (15) 

and where A9 is an angular interval to be defined. In commercially available haze-meters [isl [l6l [sij . one for this 
angular interval uses the value A9 = 2.5° 64], and for the present study this value will be adopted. Notice that the 
angles 9± can be related to the lateral momentum variable, q±, in accordance with Eq. ([6|). Furthermore, the angular 
interval A9 is related to a corresponding momentum interval in the following way 

Aq = i((7+-g_) = ye;;;" — cos e sin A6I. (16) 
2 c 

Haze, as defined above, can readily be related to the mean differential reflection or transmission coefficients {dU/d9) 

(cf. Sec. IHI Bp . In terms of these quantities, haze, that in general will depend on the angle of incidence 6*0, can be 
written in the form 

W) = (17a) 

= 1-I{9^,9+), (17b) 



I{Oa,0b) - h f^d9(^), (17c) 



where one has introduced 

^ J 9^ \ 99 

and where U has been defined earlier in Eq. In the transition from Eqs. (|17ap to (|17b|) it has been used that 

/ (— 7r/2, 7r/2) = 1 by definition. The reason for the presence of the factor \/U in Eq. (|17cp is that haze is defined in 
terms of the fraction of the reflected or transmitted light, while {dU/d9) is defined as the fraction of the incident power. 
Hence, this factor is present to ensure that the defining expression of haze has been given the correct normalization. 
From the definition, Eqs. P?|) . it follows that haze is a dimensionless number between zero and one [1^. Furthermore, 
notice that an ideal scattering system, i.e., one with no surface or bulk randomness, will correspond to a haze value 
of H{9q) = for all angles of incidence, since all light will be reflected or transmitted into the specular direction [66l |. 
However, as the randomness of the scattering system is increased, and therefore the reflected or transmitted intensities 
as a consequence become more and more diffuse, the corresponding haze value will increase. To reach a haze of one 
(Ti. = 1) is, however, rather unlikely for random systems encountered in practical situations since it will require a 
vanishing intensity over the angular interval (0_, 0+). Surface random systems with such a property can, however, be 
artificially manufactured (52| . From a practical point of view, a more likely scenario for a strongly random system is 
probably that of a Lambertian diffuser (s^, i.e., a scattering system giving raise to {dU/d9) oc cos 6* independent of 
angle of incidence. For such a scenario, haze at normal incidence will be 7i(0) = 1 — sin A6' ~ 0.956 < 1 according to 
the definition ()17p . For naturally occurring surfaces this is probably a more realistic upper limit of haze. 



B. A naive approximation to haze 



In many situations encountered in practical applications of the concept of haze, it is illuminative to have available 
an approximate expression to the formal definition (jl7p . Often the dependence on various parameters on that haze 
will depend, can be made more apparent via approximate expressions. 

Before considering an analytic approximation to haze (see next subsection) we will present some approximations 
that are based on the concept of coherent and incoherent scattering. Such picture is often useful to bear in mind 
when working with the haze (and gloss) concept. 

If the scattering system is not too diffuse, the main contribution to the /-integral present in Eq. (|17b) will come 
from the coherent component of the scattered or transmitted light. Furthermore, since the coherent component is 



non-zero outside the angular interval from 9^ to one arrives at the following approximation to haze 



n{eo) ^ 1-77 / de(—) , (ISa) 



^ J-./2 \ de / 

based exclusively on the coherent component of the scattered or transmitted field. Thus, it is natural to refer to 
Eq. (|18ap as a coherent approximation to haze. From a mathematical point of view, one may rewrite Eq. (jl8ap by 
noting that the last term is just one minus the corresponding incoherent component. By recalling Eqs. (|12p and 
one, therefore, arrives at the following equivalent incoherent approximation to haze 

1 r'"^ /dU\ 

Later it will be shown that the coherent approximation, Eq. (jl8ap . only takes into account the surface roughness, 
and not its correlation. The incoherent approximation, Eq. (|18b|) . however, will be shown to also depend on the 
surface correlation. Even though the coherent and incoherent approximations to haze are equivalent from a purely 
mathematical point of view, they give rise to different physical interpretations. Since Eq. (jl8bp includes the dependence 
on the surface correlation, as well as surface roughness, we will prefer this approximation over that of Eq. p8ap . 

An approximation is not worth much without information about its range of validity. In order to get an idea of 
the accuracy of the coherent and incoherent approximations, the Lambertian diffuser will be considered once more. 
Such a diffuser represents in many ways a worst case scenario since there is no coherent component at all in this 
particular case. Due to the vanishing coherent component, one therefore, within the approximations of Eqs. (jl8p . has 
that Ti ~1. However, above one found by using the formal definition of haze, Eqs. (fT7|). that Ti. ~ 0.956. Hence, the 
error obtained by calculating the haze from the approximate expressions (ITS)) is of the order of 4.5% for this highly 
diffusive case. For less diffusive surfaces the error is expected to be less. 



C. An analytic approach to haze 

In this subsection, an analytic expressions for haze will be derived. This will be achieved by applying the so-called 
phase perturbation theory [sj, [ssj . This approximation is reviewed in Appendix and it can, at least in reflection, 
be viewed as an extension (or correction) to the more well-known (and used) Kirchhoff approximation (23. [25. 26, 27j. 
In particular, phase-perturbation theoretical results reduce in the limit of large correlation length for the surface 
roughness to those that can be obtained from Kirchhoff theory. In addition, phase-perturbation theory naturally 
can handle transmission problems in an analytic fashion while such a generalization is not straight forward for the 
Kirchhoff approximation. 

Within phase-perturbation theory, the mean DRC or DTC, {dU/dff), can be written in the following form (cf. 
Appendix [X]): 

/ dU\ 1 £„i UJ cos^^i ^,^,2 



\ o6 / L 27rc cos 

with uo{k) being the polarization dependent Fresnel reflection or transmission coefflcients [30| corresponding to the 
scattering system of a planar (non-rough) interface, and 

/oo 
rf^e%-fe)«e-^A=(?l'c)M/(„)^ (19b) 
-OO 

where 

A(g|/fc) = ( + reflection^ 

^ ' ' 1 ai[q)~aQ[k), m transmission ^ ' 

denotes the momentum transfer perpendicular to the mean surface. In writing Eqs. (jl9p . we recall that a and W{u) 
are the surface roughness and correlation function, respectively (cf. Eqs. ([T])), and that L denotes the length of the 
rough surface measured along the mean surface. Notice, that for a planar surface (cr = 0), J{q\k) will be proportional 
to S{q — k) and hence the effect of the surface roughness in Eq. p9ap is fully contained in the J((7|fc)-integral. 

In order to derive an approximate analytic expression to haze, we start by assuming that (T'^|A'^((7|A:)| <C 1 for all 
values of q in the radiative region defined by \q\ < ^/e^uj/c. Within this approximation, the J-integral takes on the 
following form 

J{q\k) ~ Le-"'^'(«l'=) [2n6{q - k) + a^A^{q\k)g{\q - A:|)l , a^\A^{q\k)\ < 1, (20) 



where Eq. Q has been used to relate W(|a;|) to the power spectrum 5(|fc|). It should be noticed that the first 
term of Eq. (|20p is a coherent (specular) contribution, while the last term represents the lowest order contribution 
from the incoherent (diffuse) field. Moreover, the coherent term depends on the parameter that defines the surface 
roughness only through the rms roughness, a, while the incoherent contribution in addition shows a dependence on 
the correlation length a via the power spectrum (?(|A:|). 

In order to calculate haze, we need to get an expression for the (integrated) fraction of the incident power that is 
either reflected or transmitted into any angle, i.e. one is looking for an expression for U as defined by Eq. (fT3|) . The 
main effect of surface roughness is to alter the angular distribution of the light being reflected from, or transmitted 
through, a randomly rough surface. However, the amount of integrated reflected or transmitted light is much less 
sensitive to the presence of surface roughness. Notice, that this is only true if we are not close to a roughness induced 
surface resonance like e.g. those due to the excitations of surface plasmon polaritons [56l |. Another situation where 
the above assumption is known to fail is in situations where the scattered or transmitted intensity is rather low for 
the planar surface so that the presence of roughness might renormalize this result in a significant way. This is for 
instance the case for the Brewster angle phenomenon [30]. Hence, for the purpose of this study, it will be assume 
that U shows little sensitivity to surface roughness so that it can be well approximated by the planar result. With 
Eq. (|19ap and a ~ one therefore has 

U^^\uo{k)\\ (21) 

where one, as before, should use to = in refiection and to = 1 in transmission. 

By substituting Eqs. (^0]) and (PT|) into Eqs. an approximate result for the /-integrals, Eq. ()17cp . used to define 
haze, can be obtained as 

1 r 



L Jg 27r am{k) 

Jq_ 27r arn{k) 



where g± = fc ± Ag. To try to obtain an estimate for the last integral of Eq. ((22|) will now be approached. Since the 
momentum interval Ag is small (cf. Eq. ()16|) ). am{q) and K{q\k) are therefore slowly varying functions of q in the 
interval < q < q+. The power spectrum gd^l), however, can in particular for large correlation lengths, Aqa ^ 1, 
become strongly g-dependent over the interval of interest even for small values of Aq. Thus, the integrand of the last 
term of Eq. ()22p . accept the power spectrum g(|p|), can be approximated by its value at q — k (specular direction). 
Hence, one may write 



.2 a2/, i,xG'(a) Aq 



l + cr^A^(A:|fc)- 



(23) 



where 



= 7^ r dqgi\q^k\) = ^ [ ' dqg{\q\) (24) 

Aq 



2AqJg^ ^^^'^ " 2AqJ_ 



denotes the power spectrum factor calculated over a momentum interval of half-width Aq around zero momentum 
transfer q — k = 0. It is important to realize that G{a) is independent of the angle of incidence, but will, of course, 
depend on the type of power spectrum used for the surface roughness. Hence, when the type of power spectrum is 
known, G{a) can be calculated for all angles of incidence, as well as for both reflection and transmission. For an 
exponential power spectrum, defined by Eq. ([3]), the power spectrum factor becomes 

2 

G(a) = — — arctan(A(7a), (25) 
Aq 

and for a Gaussian power spectrum g{\q\) = ^y^^aex^p{—q■^a'^ /A), often used in practice, G{a) = {tt / Aq) eii{Aqa/2) 
where erf(-) is the error function [57|. Notice that whenever Aqa 1, the power spectrum factor may be expanded 
with the result that G(a) is directly proportional to the power spectrum at zero momentum; G(a) ~ 5(0). In general, 
however, this is not the case. 

In the spirit of phase-perturbation theory [s^ . [ssl . [ssj , one observes that the terms in square brackets in Eq. 
are the first few (non-trivial) terms of an exponential function. Hence, one may write 



I{e-,e+) ~ exp 



(26) 



and after substituting this expression back into Eq. (|17bp . the foUowing approximate expression for haze is finaUy 
arrived at 



1 — exp 



'a^k\k\k) 1 



G(a)Ag 



(27) 



Here the perpendicular momentum transfer, A(g|fc), has aheady been defined in Eq. (|19cp . the power spectrum factor 
G{a) in Eq. (j24p , and the momentum interval Ag in Eq. (jl6p . The approximation (|27p and its numerical confirmation 
(to be presented later), are the main results of this study. 

There are several important observations to be made from the approximate expression to haze ( I27p : First, it should 
be observed that according to Eq. (|27p haze can be written in terms of two dimensionless quantities: aA{k\k) and 
G{a)Aq. The product of the rms- roughness of the surface topography and the perpendicular momentum transfer 
of the reflection (or transmission) process, aA{k\k), does depend on the "amount" of roughness but not on how it 
is being correlated. The quantity, G{a)Aq, on the other hand, depends on the power spectrum, and is therefore 
sensitive to the type of height-height correlation function respected by the rough surface. Second, the dependence on 
the angle of incidence onZj/ enters through the momentum transfer: aA(k\k) (perpendicular) and Aq (lateral). Third, 
to apply the approximation (f?/]) in practical applications, it is the power spectrum of the surface roughness around 
zero momentum that is of interest since only this portion of the power spectrum enters the definition of the power 
spectrum factor G{a). This is important to realize since the whole power spectrum, and in p articular its tail, is often 
difficult to assess in a reliable way from direct measurements of the surface topography [26| . 

Another quantity used frequently by the optical industry to quantify visual appearance of surfaces are, as mentioned 
earlier, gloss [11]. Gloss is, crudely speaking, a measure of how specular a surface appears in reflection or transmission. 
This can be measured by e.g. the fraction of the reflected (transmitted) light that is reflected (transmitted) into in 
a small angular interval around the specular direction O. Hence, gloss is more or less complementary to haze, i.e. 
it can more or less be written as 1 — 7^(6*0). Hence, the last term of Eq. ([27]) may therefore be considered as an 
expression for gloss. In practice one talks of at least two types of gloss; wide angle gloss and specular gloss [ip]. They 
are distinguished by the values of Ad (and therefore Aq) used to deflne them. When using the expression for gloss 
described above, one has to substitute the correct values for A9 (as well as ^o)- The expression for gloss, 1 — T-L{9q), is 
probably best suited for specular gloss due to the approximation introduced in order to arrive at Eq. (|23p . However 
for large values of AO it is not expected to work too well. 

Previously, Alexander-Katz and Barrera [45^, while studying gloss (in reflection), found that the reduced variables 
for this problem were ((t/A)cos6'o and (a/A) cos 6*9. The results of these authors, have quite a few common features 
with those presented in this study. However, this is probably not so surprising since gloss and haze can be viewed as 
more or less complementary quantities. In particular, in reflection, we find that aA{k\k) indeed scales as (cr/A)cos0o 
(cf. Eq. (|29ap below). Hence, this agrees with the result found by Alexander-Katz and Barrera for gloss [45j. However, 
for the correlation function dependent quantity, we do not in general seem to agree with these authors. Only in the 
limit Aqa <C 1, when G{a) ~ 5(0)a as mentioned earlier, do we tend to get the same scaling behavior as reported 
previously by Alexander-Katz and Barrera. The reason for this "discrepancy" , is related to the level of accuracy 
applied in the approximation of the integral of Eq. ([55]). 

When |crA(A:|fc)| ^ 1, the exponential function of Eq. (I77|) can be expanded with the result that 



^(6*0) ~ cj'^K^{k\k) 1 - 



G{a)Aq 



|aA(fc|A:)| < 1. 



(28) 



This has the consequence that the ratio of haze in reflection (7^5(6*0)) and transmission (HtiOo)) will be independent 
of the rms-roughness a and only depend on the correlation length a (and the form of the power spectrum) as well as 
the parameters defining the scattering geometry. 

For completeness, the expressions for haze in terms of the "defining" quantities will be explicitly given. With 
Eqs. (HH), (fT9cll and ^ one has for reflection 



HsiOo) ~ 1 — exp —IGneoi^—j cos 
while in transmission the haze can be approximated by 

HtiOo) 0^ 1-exp 
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(29b) 



In obtaining Eq. (|29bp it has been used that for the specular direction (in transmission) cos 8t = \/ 1 ^ (^o/ei) sin 



Before closing this section, it should be stressed that the approximate expression to haze, Eq. ((?7|) and related 
expressions, are based on phase-perturbation theory. Hence, its validity will become questionable when multiple 
scattering starts to contribute significantly to the scattered or transmitted fields. We stress that phase-perturbation 
theory is not a small amplitude perturbation theory, so it may still give reliable results for strongly rough surfaces in 
the large correlation length limit {i.e. small slopes). The accuracy of Eq. will be investigated in Sec. fVl 



In addition to haze, gloss is a quantity frequently used in the optical industry to quantify optical materials. As haze 
quantifies the fraction of the reflected or transmitted light that is directed outside a given angular interval, gloss, on 
the other hand, is related to the amount of light that falls inside the same interval. 

Mathematically, gloss, G{0o), is defined as the integral I{d+,9-) of Eq. (|17cp with 9± adapted to the particular 
definition of gloss of interest. Hence, in terms of the surface parameters, gloss can be expressed by Eq. (12^ . or by 
the subsequent approximate expressions that followed it. In particular take notice of Eq. p6|) that represents an 
approximate formula for gloss 



In light of Eq. ([27|) . one, as expected, observes that 'Hi{do) + Qi{9o) = 1, where i — s (reflection), t(transmission) and 9± 
are the same for both haze and gloss. This means that when the value of haze is increasing, gloss is reduced and visa 
versa. We stress that both quantities in the equation above refer to reflection or transmission. Mathematically there 
is no problem in consider gloss in transmission, however, this term is not commonly used in practical applications. [B^I 
For the implications and validity of the above approximate expression to gloss, Eq. ([30|) . the interested reader is 
referred to a separate publication [soj . 



The approximate expression for haze, Eq. ()27p . is based on phase-perturbation theory, and is therefore not rigorous. 
In order to investigate how well this analytic approximation is performing, and when this scaling form breaks down, it 
will be compared to what can be obtained from a rigorous computer simulation approach [sli . |60| . Such an approach 
is formally exact, since it solves the Maxwell equations numerically without applying any approximations. It will thus 
take into account any higher order scattering process, and not just those accounted for by phase-perturbation theory. 

Such exact Monte Carlo simulations can be performed by formulating the Maxwell equations as a coupled set of 
integral equations. This is done by taking advantage of Green's second integral identity in the plane as well as the 
boundary conditions satisfied by the fields and their normal derivative on the randomly rough surface. These integral 
equations can be converted into matrix equations and solved for the sources — the fields and their normal derivative 
evaluated at the surface. From the knowledge of these sources, the scattered (transmitted) field at any point above 
(below) the surface may be calculated and therefrom the mean DRC (DTC). The whole detailed procedure for doing 
such Monte Carlo simulations can be found in Refs. [3l|, [gO, [Ml ■ 

In the numerical simulation results for haze to be presented below, one will first calculate {dU /d9) and then, by 
numerical integration, calculate haze directly from Eq. (jl7p . An example of a mean differential reflection coefficient 
curve, obtained by numerical simulations for a polymer material, is given in Fig. [31 In obtaining this result one 
assumed a Gaussian height distribution function of cr/A = 0.058 and an exponential correlation function of correlation 
length a = 1.58A. These parameters are similar to those found for the measured surface of Fig. [2l The vertical 
dash-dotted lines that can be seen in Fig. [3] are at an angular position 9± = 2.5° — i.e. at the angles about the 
specular direction 9s ^ 9q ^ 0° used in the definition of haze. By numerical integration one may from this result 
calculate haze, and for this particular example one finds 'Hs{9o = 0°) = 0.34. 

In Figs. [4] rigorous numerical simulation results for haze (open symbols) are presented versus various parameters 
of the scattering geometry in both reflection and transmission. Here an exponential correlation function Pr(|a:|) — 
exp(— |a;|/a), that is a reasonable fit to measured data for some polymer films (see Figs. [5]) has been assumed. 
Furthermore, the light of wavelength A = 0.6328/iTO was incident normally onto the mean surface (Figs.HJa) and (b)), 
and for the dielectric constant of the media involved, £o = 1 and £i = 2.25, were used. For all simulations, the results 
were averaged over at least = 500 surface realizations. The vertical dashed-dotted lines present in Figs.SKa) and 



D. Gloss 




(30) 



V. SIMULATION RESULTS AND DISCUSSION 



A. Comparison to Monte Carlo simulations 



(b) correspond to the roughness parameters used in obtaining the results of Fig. [3] The soHd hues in Figs. U] are the 
predictions of the Eq. ()27p — the approximate expression to haze. In particular, Fig. Hl^a) presents the dependence 
of haze vs. a/X for a fixed value of the correlation length a/A = 1.58 and the angle incidence was = 0°. This value 
for a/X corresponds to the vertical dashed-dotted line in Fig. (H^b). From Fig. [DJa) it is observed that the analytic 
approximations (solid lines) performs impressing well, also for the roughest surfaces considered. It should be noted 
that both the (cr/A)^ increase of haze for low levels of roughness (cf. Eq. (j28p ) as well as the "bend-off", or saturation, 
that takes place for the haze for strong roughness seem to be correctly predicted by Eq. (|27|) . In this latter case, 
however, the numerical simulation results saturated around Ti ~ 0.95-0.96 while the theoretical curves approach the 
value of one. Hence, for strongly rough surfaces one has a relative error of about 5% as discussed in an earlier section. 

In Fig. mb) the dependence of haze vs. correlation length a/X for a/X = 0.058 (vertical dashed-dotted line) and 
9q = 0° is depicted. The agreement, also in this cases, is rather satisfactory. However, from this figure there is an 
indication that at smaller correlation lengths, the agreement becomes less good. This is caused by phase-perturbation 
theory not being a good approximation in the small correlation length limit (that corresponds to large local slopes 
a/a). 

Finally in Figs. 111(c), how the angle of incidence influence haze is studied. Only positive angles of incidence are 
being considered since the scattering geometry is so that there is a symmetry with respect to a change of sign in Oq. 
As was done to obtain the results of Figs.HKa) and (b), one has also here fixed the roughness and correlation length to 
the values used in obtaining the results of Fig.O i.e., a/X — 0.058 and a/X — 1.58. In the case of transmission (lower 
panel of Fig. lU^c)), the agreement between the Monte Carlo simulation result and the analytic approximation is of a 
good quality for all angle of incidence considered. In reflection (upper panel of Fig. HKc)), the agreement is of a good 
quality only for the smallest scattering angles. However, as the angle of incidence approaches roughly 9o — 55°, from 
below or above, the disagreement between the simulation and approximate result (for reflection) becomes pronounced. 
The reason for this discrepancy is the so-called Brewster angle phenomenon [s^l • This phenomenon express itself for 
the planar geometry in p-polarization by the reflection coefficient being exactly zero at the Brewster angle 6b defined 
by tan^^B = ei/eo- However, as roughness is introduced into the system, the reflectivity of the (rough) surface 
will not be zero any more, not even at the Brewster angle, but will instead go through a minimum for an angle of 
incidence close to Ob (the "quasi" -Brewster angle phenomenon). This has the consequence that haze (in reflection and 
for p-polarization) will go through a corresponding maximum for the same angle of incidence. So in the region about 
the Brewster angle, the presence of surface roughness will strongly renormalize the corresponding planar geometry 
result with the consequence that the integrated reflected energy is not any more well approximated by Eq. (|2T|) . By 
taking into account roughness in the estimation of the total integrated scattered intensity, one will most likely be able 
to also predict the behavior of haze for such angles of incidence with more confidence. However, the penalty for doing 
so, is that the resulting expressions become much more complicated and must be evaluated numerically. Since this is 
not the aim of the present study, we will not follow up this line of actions here, but only keep in mind that the simple 
approximation (j27p breaks down around the Brewster angle. Notice that there is no Brewster angle phenomenon in 
transmission as shown explicitly in the lower panel of Fig. HJc), nor is there any such phenomenon in reflection (or 
transmission) for s-polarized incident light. Hence, for s-polarization, the approximate expression for haze, Eq. (j27p . 
should apply for all angles of incidence, something that has been confirmed by numerical simulations (results not 
shown) . 

In Figs.[3Ka) and (b), one, or more, of the parameters that characterize the surface roughness were fixed to constant 
values. It is, however, important to get a more complete picture of the quality of our analytic expression to haze. 
This can be achieved by allowing both a/X and a/X to vary freely (within certain limits). Figs. [5] depict contour 
plots for the variations in haze vs. both of the two above mentioned parameters in reflection and transmission. 
The figures in the left column. Figs. [IJa) and (c), show contour plots of haze as obtained by rigorous Monte Carlo 
simulations. In the right column, i.e. in Figs. E^b) and (d), the corresponding plots obtained from the analytic haze 
(approximate) expression, Eq. (j27p . are presented. By in Figs.[S]comparing the numerical simulation results to those of 
the corresponding analytic predictions, one can conclude that the quality of the analytic expression (|27p is remarkably 
good over large regions of parameter space. This is particularly the situation when considering reflection. As a general 
trend, it seems fair to say that the approximation (j27|) performs the best for large correlation lengths and smallest 
rms-roughness. This is in particular the case for haze in transmission for which the approximation is poorer than in 
reflection. These findings fit the picture that phase-perturbation theory can be looked upon as a generalization of the 
more familiar Kirchhoff approximation that is known to work the best in the large a/X limit [l^jISo]- Moreover, notice 
that when the large correlation length limit is taken, for fixed rms-roughness, the small slope limit is approached since 
the average slope of the surface is proportional to a /a. 

So far, a semi-infinite dielectric transparent medium bounded to vacuum by a randomly rough interface has been 
considered. From a practical point of view, it will be rather interesting to also investigate how haze depends on the 
surface parameters for, say, a film geometry. For such a case, the derivation of an approximate expression for haze, 
analogous to those presented in Sec. IIVI for a semi-infinite medium, is lengthy, but can be performed. However, such 



expressions will not be presented here. Instead we would like to add that rigorous numerical simulations for a film 
geometry have been preformed for normal incidence. The results of such simulations show that the characteristic 
dependence of haze on the parameters cr/A and a/ A, originally found for a semi- infinite randomly rough surface also 
seems to hold true for the film geometry. 

B. Comparison to experimental results 

So far in this paper we have mainly considered one-dimensional surfaces. However, naturally occurring surfaces, 
as well as man-made surfaces generated by, say, an industrial process, are usually two-dimensional. Neither is it not 
uncommon that their statistical properties are not well described by simple mathematical distribution and correlation 
functions of the form often assumed in theoretical studies. It is therefore an open question if the approximate result for 
haze, Eq. (|27p . obtained for one single randomly rough one-dimensional surface has any relevance for the more com- 
plicated scattering systems encounter in practical applications. In the most general case our approximate expression 
is obviously not suitable. One may still hope, however, that the general behavior of haze found in Eq. ((27|) could be 
taken over to higher-dimensional and more complicated scattering geometries. In particular, for weakly rough isotropic 
surfaces at normal incidence there are hopes that a one-dimensional approach might work reasonably well. This is 
so since under such circumstances the (two-dimensional) mean differential reflection and transmission coefficients 
are rotational symmetric, i.e. no ^-dependence. Consequently, a one-dimensional mean differential or transmission 
coefficient might be enough to catch the main angular dependence of the scattering up to cross-polarization effects. 
The purpose of this sub-section is to look into these questions, and to compare experimental measurements with what 
can be obtained from a one dimensional computer simulation approach of the type applied previously in this paper. 

To investigate this further, we will consider a melt blown film of linear low density polyethylene (LLDPE). The 
film and LLDPE type was referred to as narrow molecular weight distribution metallocene, and described in detail, 
in Ref. [4]. The haze of the film was measured by a spherical haze-meter (Diffusion System, type M57) according to 
the standard [isj . and (in transmission, at normal incidence) one obtained for this material a haze of 18.9%. Light 
scattering at the surface, as well as in the bulk, contributed to this haze value. The bulk contribution was estimated 
by measuring the haze of films coated by glycerol on both sides. Glycerol has a refractive index of 1.47, which matches 
closely the refractive index of the blown films (approximately 1.5). Under this condition the haze was significantly 
reduced (to 5% for the film mentioned above, leaving, crudely speaking, about 14% haze resulting from the surface 
roughness). The films that were experimentally studied in Ref. 4] (with that of Fig. [2] being among them) gave a 
haze in the range 12-30%. When embedded in glycerol, the corresponding values were reduced to 3-6%. The film 
with the highest haze was measured to be 31.9%, and reduced to about 1/10 of this value when submerged in glycerol 
(3%). This indicates that the contribution to haze from surface roughness dominates over the contribution from the 
bulk. This is consistent with the assumption made in the theory section of this paper (cf. Sec. IIII[) . 

The surface topography of the above mentioned polymer film was measured and served as the basis for the char- 
acterization put forward in Fig. [D The solid lines in Figs. [11(b) and (c) represent a Gaussian (Fig. [11(b)) and an 
exponential fit (Fig[2l(c)) to the height distribution and height-height correlation function respectively. These func- 
tions are characterized by a root-mean-square roughness of cr = 0.04 and a correlation length of a = 1.3 /im. 
Notice that the functional fits performed in Figs. [1Kb) and (c) are of reasonable quality and in particular for the 
height distribution function. In the one-dimensional Monte Carlo simulation results to be presented below, one has 
used the functions represented by the solid lines in Figs.[lKb) and (c) as a basis for generating the underlying ensemble 
of surface realizations. 

In Fig. [HI the logarithm of the experimentally obtained angular intensity distribution, log 7(6', 0), of light being 
transmitted through the polymer film system described above is depicted. In obtaining these results, a HeNe-laser at 
wavelength A — 0.6328/im was used as a source for the unpolarized, normal incident light (6*0 — 0°). As expected, a 
strong specular transmission peak as well as a large dynamical range in intensity (almost seven orders of magnitude) 
are observed. Moreover, a weak anisotropy in the angular distribution of the transmitted light can be observed in 
Fig. [H] as represented by the horizontal line of enhanced intensity. This anisotropy is caused by the polymers being 
partially oriented along the direction of the flow in the production of the polymer film. In Fig. [6] this direction 
corresponds to the vertical. Except from the weak anisotropy the angular distribution of the transmitted light is 
rather isotropic. 

In Fig. [3 a horizontal cut through the center of Fig. [H] is represented by open symbols. The solid line in this same 
figure represents the rigorous (one-dimensional) Monte Carlo simulation results for the mean differential transmission 
coefficient. In obtaining this latter result a film geometry of mean thickness d = 40/iTO was used where the uncorrelated 
upper and lower rough (one-dimensional) interfaces were described by the parameters that were derived from the 
measured surface topography of the film (see Figs.[l|). In particular, these surfaces were characterize by the functions 
represented by the solid lines of Figs. [Hb) and (c) (see also the caption of Fig. [7| for the surface parameters). 



From Fig. [7] it is observed that there is a quite reasonable agreement between the measured and simulated angular 
distribution of the transmitted light. It is at the largest angles of transmission that the discrepancy starts to emerge, 
while the angular distribution in the central part, that represents the main part of the transmitted energy, seems to 
be well accounted for by the simulation result. For the largest transmission angles the Monte Carlo result seems to 
underestimate the transmitted power. This situation is in fact not unexpected; In the experimental measurements 
scattering from the bulk is also present, while such effects has not been taken into consideration in the approach used 
to produce the solid line of Fig. [71 It is, in fact, well known that bulk scattering tends to enhance the transmitted 
power into large angles of transmission |62| . 

As mentioned above, one by direct measurements found the haze in transmission for normal incidence to be 
'^215(0°) ~ 0.189 for the experimental sample (with both surface and bulk randomness). From the Monte Carlo 
simulation results one on the other had found a haze of 7iii3(0°) ~ 0.054. However, in order to be able to compare 
this result with the experimentally available value (7i2_D(0°)), one has to introduce a (multiplicative) constant that 
accounts for the difference between the one- and two-dimensional geometry (assuming normal incident light). Doing 
so results in a haze of about 14-15%, a result that compares favorably to the experimental haze value resulting from 
surface randomness. Hence, one has demonstrated that the one dimensional Monte Carlo approach can be applied to 
predict reasonably well the haze of an isotropic experimental sample at normal incidence. This is indeed encouraging 
results, but further work is needed, however, to determine the full region of applicability of the Monte Carlo approach. 

VI. CONCLUSIONS 

In conclusion, we have studied the dependence of haze and gloss with the parameters that normally are used to 
characterize randomly rough surfaces — the rms-roughness a and the height-height correlation length a. Based on 
phase-perturbation theory, we have derived analytic expressions that represent approximations to haze and gloss 
for a one-dimensional Gaussian rough surface. It is demonstrated that haze(gloss) increases(decreases) with cr/A 
as exp(— A((t/A)^) and decreases(increases) with a/A in a way that depends on the specific form of the correlation 
function being considered. This latter dependence enters into the expression for haze and gloss as the average of the 
power spectrum taken over a small interval, Aq, around zero momentum transfer. In the limit Aqa <^ 1 one obtains 
that the haze and gloss depend linearly on the correlation length. 

The range of validity for these approximate expressions to haze and gloss put forward in this paper were assessed 
by rigorous numerical Monte Carlo simulations. They were found to agree remarkably well over large regions of 
parameter space. Furthermore, some experimental results for the angular distribution of the light being transmitted 
through a polymer film was presented. It was found to fit reasonably well with the prediction from the Monte Carlo 
simulations, and consequently the predicted value of haze was found to agree relatively well. 
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APPENDIX A: PHASE PERTURBATION THEORY 

In this appendix we will describe the so-called phase perturbation theory, as well as deriving some analytical ex- 
pressions, on which the main text relays. What has become known as the phase-perturbation theory, was originally 
developed by J. Shen and A. A. Maradudin [s^] for non-penetrable media. The method was later extended to one- 
dimensional randomly rou gh p enetrable media by Sanchez-Gil et al. [s^. The method has also been formulated in 
an explicit reciprocal way [55|. In phase-perturbation theory it is the phase of the field that is determined pertur- 
batively [1^, and it has proven well suited for reflectivity studies [s^, [63|. One of the interesting features of this 
pertubativc method is that in the large limit of a/ A, with a being the surface correlation length and A the wavelength 
of the incident light, it reduces to the more well-known Kirchhoff approximation [l^j H^j [26,, J2J ^^r a non-penetrable 
medium. Furthermore, as a becomes comparable to A, phase-perturbation theory represents a correction to the Kirch- 
hoff result. An additional practical advantage of phase-perturbation theory is that it remains its analytic form also 
for penetrable or absorbing media, something that is not the case for the Kirchhoff approximation. 

Let us start our discussion of phase-perturbation theory and what can be derived from it, by letting U{q\k) col- 
lectively denote either the reflection or transmission amplitudes, Ry{q\k) or T^{q\k), introduced in Sec. IIII Bl By 



taking into account the boundary condition at the rough surface, and assuming that there is no down-going (up- 
going) scattered (transmitted) waves even close to the rough interface [68] a single integral equation for U{q\k) can 
be derived [H, [3l| (the reduced Rayleigh equation). Based on this integral equation one can derive the following 
expression for amplitude 

U{q\k) = uo{k) / e-^(«-'=)^e^^(«l'=)'^(^) (Ala) 



where uo{k) is the Fresnel coefficient for the corresponding planar geometry. In writing this expression we have 
introduced 

A(g|fc) = {""['^J+^^SS' . (Alb) 

^ ' ' 1 ai(q) — ao(k), m transmission ^ ' 

It should be noticed that if we were dealing with a more complicated scattering geometry then the one depicted in 
Fig. [U say a film geometry, the amplitude could still be expression in the form (jAlaP if only one of the interfaces were 
randomly rough. In this case A(g|fc) will take on another form then the one given above. 

In order to make contact with observable we will need an expression for |J7(g|A:)p. In fact, more precisely it is the 
average over this quantity that we should be interested in since the surface is randomly rough. With Eq. (|Ala|) . we 
find 

{\U{q\k f) = |^.o(fc)|' r dx r dx' e-'(«-'=)(---') /e»A(,|fe)[c(.)-C(.')]\ ^ (A2) 

where the average, denoted by (•••), is assumed to be taken over an ensemble of surface realizations of the surface 
profile function C,{x). Furthermore, we have here assumed that K(ci\k) is real (or close to being real), as it will be 
automatically in reflection. 

With the change of variable x' — x + u the above equation becomes 

(\U{q\k)\^) = \uo{k)\^ ( dx f du 6^9-'=)" /e*^(«l'=)^^("A , (A3) 

with AC(u) = (,(x) — C,(x^u). If C(x) is a stationary random process, then the average in Eq. (|A3p will be independent 
of x. Therefore, the x-integration in the same equation will give the contribution L, that is the length along the x- 
direction of the rough surface. Furthermore, if in addition to being stationary C also is a Gaussian random process, 
as we assume here according to Sec. [TTl then A({u) will also be a Gaussian random variable. Hence, the average in 
Eq. (|A3|) can be done analytically with the result that 

{\U{q\k)\') = L\uo{k)\^ d^e*(«-fc)«e--'A^(9l'=)[i-^Hl 

= \uo{k)\^ Jiq\k), (A4a) 

where 

J{q\k) = Le-"'^'*^!'^) / d,,e*(9-'=)"e"'^'(«l^-)^("), (A4b) 



with W{u) being the surface correlation function as defined in Sec. [TTl 

With Eqs. fXt)) and (jllbp . one for the mean DRC or DTC, collectively denoted {dU/d9), finally obtains 

/ — \ = i^^^^l.o(fc)|^J(9|fc). (A5) 
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FIG. 1: The scattering geometry used in tliis study. The rough surface is defined hy z = Ci^)- The region above the surface, 
2 > C(a;), is assumed to be vacuum (eo(t^) = 1), while the medium below is a dielectric characterized by a frequency-dependent 
dielectric function ei(a;). Notice for which direction the angle of incident {6o), scattering (Os), and transmission {6t) are defined 
as being positive. An angle of transmission is only well-defined if the lower medium is transparent, i.e. if i?eei(a;) > 0. 
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FIG. 2; (Color online) Experimental height measurements obtained from Atomic Force Microscopy (AFM) measurements on 
a linear low density polyethylene (LLDPE) film. The material was of the Narrow (molecular weight distribution) Metallocene 
grade, (a) A contour plot of the experimental height function measured over a 50 x 50 {/J-rn)^ quadratic area. The color code is 
so that red corresponds to 0.015 fj,m and blue to —0.015 /xm. (b) The height distribution function P(C) calculated from the AFM 
data (open circles) and a Gaussian fit (solid line) corresponding to cr = 0.04 /im. (c) The height-height correlation function 
W^(|a;|) obtained from the AFM data (open circles) and fitted with an exponential correlation function W(|a;|) = exp(— [a;|/a) 
(solid line) corresponding to a correlation length of a = 1.3 /im. 
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FIG. 3: (Color online) The mean differential reflection coefficient, (dRp/dOs) , vs. the scattering angle 6s as obtained by rigorous 
Monte Carlo simulations. The incident p-polarized light of wavelength A — 0.6328/im was incident at normal incident {0o =0°) 
onto the rough dielectric surface, and for the incident wave a finite sized beam of half-width g = 6.4 was used in order to 
reduce end effects. This surface of length L = 25.6fim — 40. 5A separates vacuum, above the surface, from the dielectric medium 
of dielectric constant £1(01) = 2.25 below the surface. The statistical properties of the randomly rough surface was characterized 
by a Gaussian height distribution function of (rms) width a = 0.037/im — 0.058A and an exponential correlation function of 
correlation length a = l^m = 1.58A. The surface was discretized at A'^^ — 500 equally distributed points, and the result was 
averaged over A^<; = 5000 surface realizations. The vertical dash-dotted lines are at 0± — ±2.5°. Notice the specular (coherent) 
peak around 6s = ^o- The haze in reflection for this surface is according to the numerical simulations Ti.{6o) = 0.33, while the 
prediction of Eq. ^ is n{do) = 0.34. 
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FIG. 4: (Color online) Haze, Ti, as a function of (a) the surface roughness a/A for a/A = 1.58 at normal incidence, (b) the 
correlation length a/A for a/A = 0.058 at normal incidence, and (c) the angle of incidence 6o for a/A = 0.058 and a/A = 1.58. 
For all figures the wavelength of the p-polarized incident light was A = 0.6328 ^m. The open symbols are results of rigorous 
Monte Carlo simulations, while the solid lines are the predictions of Eq. (|27|l . The dashed-dotted line in the upper panel 
(reflection) of Fig. |4jc) corresponds to the position of the Brewster angle, Oo = 6b determined by tan^ 6b ~ ei/eo- 



(a)Rigorous simulation result in reflection 



(b)Analytic approximation in reflection 
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(c) Rigorous simulation result in transmission 



(d)Analytic approximation in transmission 
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FIG. 5; (Color online) Contour plots of the rigorous Monte Carlo simulation results (Figs. [5ja) and (c)) and the analytic 
approximation (|27p (Figs.[Sl[b) and (d)) for haze obtained in reflection and transmission for light of wavelength A = 0.6328/im 
incident normally (9o — 0°) onto the rough surface. The dielectric media that was separated from vacuum by a rough interface, 
was characterized by the dielectric constant ei — 2.25. AH results were averaged over at least N^^ — 500 surface realizations. The 
random surfaces were all characterized by a Gaussian height distribution function of standard deviation a, and an exponential 
height-height correlation function of correlation length a. Overall the agreement between the analytic and rigorous simulation 
results is satisfactory over large regions of parameter space. 



FIG. 6: (Color online) Experimental results for the angular distribution of the light being transmitted through a melt blown 
LLDPE film of the type described and partly characterized in Fig. [S] The (mean) thickness of the film was d = 50fim, and 
the wavelength of the light being incident normally onto the top mean surface was A = 0.6328/im. It is the logarithm of the 
transmitted intensity in arbitrary units that is presented. The measurements were conducted with a spectro-photo-goniometer 
built by SINTEF. The weak anisotropy seen in the transmitted intensity is caused by the polymers being preferentially oriented 
in the flow direction (vertical direction in the figure). 
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FIG. 7: (Color online) The mean differential transmission coefficient (dTp/dOt) vs. angle of transmission 9t for light of 
wavelength A — 0.6328 ^m impinging normally { 80 = 0°) onto a d = 40^m thick LLDPE film (ei = 2.25). The experimental 
results (open circles) corresponds to a cut through Fig.jH] Due to the use of arbitrary units in the experiment, the amplitude of 
the measurements was adjusted to fit that of the simulation results. In obtaining the simulation results (solid line), the surface 
parameters represented by the solid lines in Figs. [2] were used, i.e. the parameters used were a = 0.04/im and a = 1.3/im for 
the exponentially correlated rough surface. In order to replicate the unpolarized incident light used in the experiment, the 
simulation results were averaged over the s- and p-polarized results. The thickness of the film was also here d = 40^m. 



